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COMPETING PARTICLE SYSTEMS EVOLVING BY 
INTERACTING LEVY PROCESSES 1 

By Mykhaylo Shkolnikov 

Stanford University 

We consider finite and infinite systems of particles on the real 
line and half-line evolving in continuous time. Hereby, the particles 
are driven by i.i.d. Levy processes endowed with rank-dependent drift 
and diffusion coefficients. In the finite systems we show that the pro- 
cesses of gaps in the respective particle configurations possess unique 
invariant distributions and prove the convergence of the gap processes 
to the latter in the total variation distance, assuming a bound on the 
jumps of the Levy processes. In the infinite case we show that the gap 
process of the particle system on the half-line is tight for appropriate 
initial conditions and same drift and diffusion coefficients for all par- 
ticles. Applications of such processes include the modeling of capital 
distributions among the ranked participants in a financial market, 
the stability of certain stochastic queueing and storage networks and 
the study of the Sherrington-Kirkpatrick model of spin glasses. 

1. Introduction. Recently, invariant distributions for the gaps in a par- 
ticle system on the real line have received much attention. In the continuous 
time setting such questions are motivated by the study of gaps in ordered 
Brownian particle systems with rank-dependent drifts and diffusion coef- 
ficients. The latter arise in the modeling of the capital distribution in a 
financial market (see, e.g., [3, 8, 11] and [14]) and as heavy traffic approxi- 
mations of queueing networks (see Section 5 of [14] for the correspondence 
between the processes of gaps and reflected Brownian motions and, e.g., 
[12, 13] and [27] for the heavy traffic approximation of queueing networks 
by the latter). In the discrete time setting these questions appeared in the 
study of the Sherrington-Kirkpatrick model of spin glasses and were allowed 
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to characterize the quasi-stationary states in the free energy model (see [22] 
for this result and [1, 24] for its extensions). In these papers the increments 
added to the particles in each time step were either assumed to be i.i.d. or 
normally distributed. The latter evolutions by Gaussian increments can be 
viewed as a time discretized version of the Brownian particle systems. 

The present paper studies the processes of gaps in finite and infinite par- 
ticle systems on the real line and half-line evolving in continuous time and in 
which particles are driven by i.i.d. Levy processes with jumps. In the finite 
case the latter are endowed with rank-dependent drift and diffusion coeffi- 
cients. We will refer to these evolutions as competing particle systems. More 
precisely, the competing particles are indexed by i G / with I = {1, . . . ,N} 
in the case of finitely many particles and I = N in the case of infinitely many 
particles. The evolution on the whole line will be referred to as the unreg- 
ulated evolution and the evolution on the half-line [b, oo) for some barrier 
b G M. as the regulated evolution. The particle systems are defined as the 
unique weak solutions (see Section 2.1 and Proposition 3.1 for more details) 
to the following stochastic differential equations: 

dXi(t) = l{Xi(t)=x {j) (t)}8j dt + Y 1 {x % (t)=x (j) (t)}^j dBi{t) 
jei jei 

(1.1) 

+ dLi{t) + dRi{t) 
for all i G I where Ri(t) = 0, i G I in the unregulated case and 
(1.2) R l (t)= £ l {Xi(s _ )+Aii(s)<b} (6-X i ( S -)-AL l ( s )) + A (iib) (t), 

0<s<t 

i G / in the regulated case. Hereby, Xm (t) < Xr 2 ) (<)<■•• is the ordered vec- 
tor of particles in which ties are broken according to an arbitrarily specified 
ordering of the initial configuration, B(t) = (Bi(t),i £ I) is a collection of 
i.i.d. standard Brownian motions, L(t) = (Lj(i),i G /) is an independent col- 
lection of i.i.d. pure jump Levy processes each making finitely many jumps 
on any finite time interval. For each i G / and s > the term ALj(s) de- 
notes the jump Li(s) — Lj(s— ) of the process Lj at time s and each Any\(t) 
is the local time process of Xi at b. To avoid any confusion we choose the 
same normalization of the local time processes as in [3, 4] and [14] so that 
^-(i,b) (t) i s one °f the summands in the semi-martingale decomposition of the 
process i|Xj(t) — b\ rather than \X{(t) — b\. Moreover, we impose the initial 
conditions Xi(0) < X 2 (0) < ■ ■ ■ and = R^O) = R 2 (0) ■ ■ ■ and, in addition, 
X\ (0) > b in the regulated case. To distinguish between the two evolutions 
we denote the particle configuration at a time t > by X(t) = (Xi(t) :i& I) 
in the unregulated particle system and by X R (t) = (X R (t):i G /) in the 
regulated particle system. 
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Heuristically, particle i in the unregulated particle system is driven by 
a Levy process whose drift and diffusion coefficients change according to 
its rank in the particle system. The regulated particle system is defined 
similarly, except that whenever a particle hits the barrier b or jumps over b, 
it is reflected in the former case and it is placed at b in the latter case. 

To simplify the notation we set Yi(t) = X^(t) and Y R (t) = X R Jt), re- 
spectively, for alH € /. Next, we can define the processes of gaps, which will 
be the main quantities of interest, by 

(1.3) Z(t) = (Z 1 (t),Z 2 (t), ...) = (Y 2 (t) - Y 1 (t),Y 3 (t) - Y^t), . . .) 
and 

(1.4) Z R (t) = (Z R (t), Z R (t\ ...) = (Y R (t) - b, Y 2 R (t) - b, . . .), 

respectively. For / = {1, . . . , N} we can set Ci(t) = e Xi ^ and note that the 
properties of the gap process Z(t) correspond to the properties of the mass 
partition 

%)(*) 1 



(1.5) 



C (1) (t) + • • • + C {N) (t) £f =1 e*0) (')-*«) (*) ' 



1 < i < N and the analogous statement is true for the process Z R (t). 

In models for the capital distribution in a financial market the processes 
Ci(t), . . . , Cjv(i) stand for the capitalizations of the different firms and the 
particles represent the logarithmic capitalizations of the firms. Due to this 
interpretation, we will call the diffusion coefficients from now on volatilities 
or volatility coefficients of the particles. In this context it is natural to al- 
low the particles to have jumps corresponding to such events as mergers of 
firms, shares emissions or jumps in the stock price. Indeed, since Merton's 
seminal paper [17] and a vast number of works that followed, it is widely 
agreed that market models allowing for discontinuities in the stock price are 
able to better account for sudden large changes in the values of the stocks, 
the incompleteness of financial markets and the imperfectness of the hedg- 
ing strategies for options when compared to models assuming continuous 
trajectories of the stock prices. Hereby, models in which the jumps originate 
from pure jump Levy processes play a central role. For a recent summary 
of the research in this direction and a list of references we refer to the book 
[9]. 

The models on the half-line describe markets in which either firms are 
bailed out if their capitalization falls below a certain predetermined barrier 
or in which capitalizations are constrained not to overcome a certain barrier, 
depending on the interpretation of the model. In the absence of jumps such 
models can be analyzed by exploiting their relation with semi-martingale 
reflected Brownian motions in the sense of [20, 25] as was done in [14] for 
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the models on the whole line [see the proof of Lemma 2.1(b) below for the 
connection]. In [16] the authors considered a related model for a market 
with two firms defaulting if their capitalization crosses a fixed threshold. 
They were able to determine how to split the total drift, representing the 
total amount of tax cuts or subsidies, between the two firms to maximize 
the probability that both of them survive. 

The mass partitions above correspond in this context to the market 
weights of the ranked market participants. For this reason, all results on 
the gap processes in the finite particle systems translate directly into cor- 
responding results on the processes of market weights. In particular, the 
invariant distributions, the existence and uniqueness of which we prove for 
the processes of gaps in the finite systems both on the line and on the half- 
line, stand for capital distributions among the ranked market participants 
which are not changing under the evolution. We also show that the gap pro- 
cesses converge to the respective unique invariant distributions in the total 
variation distance. This corresponds to the statement that the considered 
financial market approaches a stable capital allocation in a strong sense. In 
addition, under the invariant distributions for the gap processes, all gaps 
are finite by definition which implies that the limiting stable capital distri- 
butions are nondegenerate, in the sense that the market weights of all firms 
are positive almost surely. 

In the context of stochastic queueing and storage networks, our model 
is closely related to the so-called Levy networks which are described math- 
ematically by Levy processes regulated to stay in an orthant by normal 
projections onto the boundary and normal reflections at the boundary (see 
Chapter IX of [2] and the references therein for more details). In particular, 
the techniques we use to prove that the processes of gaps possess unique 
invariant distributions and converge to the latter can be used without al- 
teration to show the same statements for an integrable Levy process with 
a nondegenerate Brownian part, whose components are endowed with neg- 
ative drifts and which is regulated to stay in an orthant, provided that its 
jumps are dominated by the drifts in an appropriate sense. The convergence 
in total variation of the process to its invariant distribution proves the sta- 
bility of the corresponding queueing network, meaning that the joint law of 
the workload processes converges in the strong sense. 

Throughout the paper we make the following assumptions. 

Assumption 1.1. (a) In all cases we assume that the volatilities a\ , oi , . . . 
are positive and it holds that 
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for some M > 1 in the unregulated case and, in addition, > Si = 82 = ■ ■ ■ 
in the regulated case. 



Our results for the systems with finitely many particles are summarized 
in the following two theorems. 



Theorem 1.2 (Invariant distributions, finite systems), (a) If I = {1, . . . , N} 

and Li satisfies 



(1. 



E 



E \ AL 



0<s<l 



< — • min(5i - S 2 , . . ■ , S N -i - S N ), 



then the gap process Z of the unregulated particle system has a unique in- 
variant distribution. 

(b) If 1 = {1, . . . ,N} and Li satisfies 



(1.9) E 



E \ AL ^ 



0<s<l 



< — -mm(-Si,Si -S 2 ,...,S N -i -S N ), 



then the gap process Z R of the regulated particle system has a unique invari- 
ant distribution. 



Theorem 1.3 (Convergence, finite systems), (a) If I = {1, . . . ,N} and 
condition (1.8) holds, then the process Z(t), t>0 converges in total variation 
to its unique invariant distribution. 

(b) If I = {1, . . . ,N} and condition (1.9) holds, then the same statement 
is true for the process Z R (t), t > 0. 



We remark that in the absence of jumps, condition (1.8) simplifies to 
Si > ■ ■ ■ > Sn. In particular, the latter is stronger than condition (3.2) of 
[14] (used there to establish the existence and uniqueness of an invariant 
distribution for the gap process of the finite unregulated system without 
jumps) which in our notation reads A X/j=i ~ Ji 

This shows that one cannot expect the conditions in Theorems 1.2 and 1.3 
to be sharp in general. 

For the infinite regulated system we show the following proposition. 



Proposition 1.4 (Tightness of the infinite regulated system). /// 
and 



(1.10) 



E 



E \ AL ^ 



< -81, 



L 0<s<l 
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then the family Z (t), t>0 is tight on Wl equipped with the product topology 
for any initial condition in 



The proof of Proposition 1.4 relies heavily on the fact that under As- 
sumption 1.1, the evolution of each particle in the infinite regulated system 
is given by a Levy process with negative drift regulated to stay on [b, oo). For 
this reason it does not carry over to the setting of the infinite unregulated 
system. At this point we also note that in the special case of the infinite reg- 
ulated system with equal drifts 61 = 82 = ••• , equal volatilities o~\ = a<i = ■ ■ ■ 
and no jumps, the processes X^(t) — b are independent reflected Brownian 
motions on R + . Due to the negativity of the drifts this implies that for each 
i £ I the process X^(t) — b converges in law to an exponential random vari- 
able. Moreover, a sequence of i.i.d. exponential random variables is almost 
surely dense in R + by the second Borel-Cantelli lemma. This shows that, 
in general, one cannot expect the infinite regulated system to have an in- 
variant distribution, since already in the described special case the number 
of particles on each nonempty interval of the form [b, y) will tend to infinity 
for any initial particle configuration. In contrast to this, in [22] the authors 
show the existence of an infinite family of invariant distributions for the gap 
process in the corresponding particle system on the whole line. In the more 
general case of nonconstant drifts or volatilities the questions of existence 
and uniqueness of invariant distributions for the gap process of the infinite 
regulated system and of the convergence of the latter to an invariant distri- 
bution are open. The approach used in the analysis of the finite systems may 
apply to this case as well. However, one has to establish the recurrence and 
irreducibility structure of the gap process in an infinite-dimensional space 
where the analysis of the corresponding reflected Brownian motion is more 
intricate. 

In [19] the authors treat the infinite unregulated system without jumps 
with the drift sequence <5i,0, 0, ... and a constant and positive volatility 
sequence o~\ = a 2 = ■ ■ ■ and are able to find an invariant distribution for the 
gap process on a subset of W. In [22] the authors characterize all invariant 
distributions for the process of gaps in the case of constant drift and volatility 
sequences. As in the regulated case the questions of existence and uniqueness 
of invariant distributions and of the convergence of the gap process to the 
latter for a general drift or volatility sequence are open. The main obstacle 
hereby is the same as in the case of the infinite regulated system (see the 
end of the previous paragraph) . 

It is a natural question to ask to provide an explicit description of the 
invariant distributions in Theorem 1.2 in the presence of jumps. However, 
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even in the case of the regulated particle system with a single particle, where 
the gap process is a one-dimensional Levy process regulated to stay nonneg- 
ative, the invariant distribution is not known explicitly. To the best of our 
knowledge the most explicit result in this direction is obtained in [5], Section 
3. There, in the described special case with the additional assumption that 
the magnitude of the negative jumps of the driving Levy process is exponen- 
tially distributed, an expression for the Fourier transform of the invariant 
distribution is given in terms of the characteristic exponent of the driving 
Levy process and the expected accumulated local time of the gap process 
at zero over a unit time interval when the initial value of the gap process 
is chosen according to the unique invariant distribution (see [5], equation 
(3.11)). A completely explicit description of the invariant distribution for 
the gap process, even for specific examples of a regulated or an unregulated 
system with multiple particles in which jumps are present, seems to be out 
of reach at the moment. 

Previous papers on Brownian systems with rank-dependent drifts and 
volatilities use extensively the results of [26] on reflected Brownian motions 
and classical results on constrained diffusion processes to obtain the exis- 
tence and uniqueness of the invariant distribution and the convergence of 
the gap process to the latter. Due to the presence of jumps, these tools do 
not apply here and are replaced by more general techniques from the ergodic 
theory of Harris recurrent Markov processes, as in [18]. More precisely, to 
obtain the results of Theorems 1.2 and 1.3, we first find the dynamics of 
the gap processes in Lemma 2.1 below. We use it subsequently to prove the 
tightness of the processes of gaps by dominating each of them by an appro- 
priate Harris recurrent Markov process. Next, we show the Feller property 
for the gap processes to establish the existence of the invariant distributions 
for the latter. Finally, to prove the uniqueness of the invariant distributions 
and the convergence of the processes of gaps, we employ the general theory 
of Harris recurrent Markov processes. 

The paper is structured as follows. In Section 2 we treat the particle sys- 
tems with finitely many particles. We show first their existence and explain 
some of their properties in Section 2.1 and then prove Theorems 1.2 and 
1.3 in Section 2.2. Section 3 deals with systems of infinitely many particles. 
We show that the latter exist and are well defined for appropriate initial 
conditions in Section 3.1 and prove Proposition 1.4 in Section 3.2. 

2. Systems of finitely many particles. 

2.1. Existence and properties of the processes. Throughout this section 
we deal with the two evolutions of finitely many particles, that is, we set 
I = {1, . . . , N}. The existence and uniqueness of a weak solution to the un- 
regulated version of (1.1) can then be seen as follows. As remarked in Section 
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2 of [3], the main result of [6] implies that for any x £ M N with x\ < • • • < xn 
there exists a unique weak solution to 

N N 

(2.1) dX-' x (t) = l {x c,* (t)=x c,x (t)} <5j + 1 {X i 0,a (t)=xg*(t)}°j dBi{t), 

(2.2) X c ' x {0)=x 

defined on some probability space (Q x , J rx ,F x ). Next, let (Q x ' 1 , J^' 1 ,W X > 1 ) , 
I € N, be copies of (fl x , T x ,¥ x ) such that on each of them a process X c,X)l of 
the same law as X c,x is defined. Moreover, let (fi L , J- L ,¥ L ) be the probability 
space on which the processes L\{t), . . . ,L^{t) are defined. Then a weak 
solution to (1.1) can be defined on the product space 

(2.3) (0,.F,P) = (jjn x ' 1 x n £ l ^7 a ' i 8J i ,(g)P s ' I «iP , 'j. 

Denoting by T\ < < ■ ■ ■ the jump times of the M^-valued process L{t) = 
(Li(t), . . . , Ljv(i)) we se ^ 

(2.4) X(0) = (X l (0),...,X N (0)), 

(2.5) X(t)=X c ^ T ^ k (t-T k ), T k <t<T k+1 ,k>0, 

(2.6) X(T fc ) = X(T fc -) + AL(T fc ), k>l, 

where we have used the notation To = 0. This gives a weak solution to the 
unregulated version of (1.1) by construction. Noting that for any other weak 
solution the joint distribution of its continuous part, its jump times and its 
jump sizes have to coincide with the corresponding quantity of the solution 
just constructed, we conclude that the weak solution is unique. Here we have 
used the uniqueness of the weak solution to (2.1), (2.2). 

The regulated process X R can be constructed on the same probability 
space as X as follows. We start with the desired initial value X R (0) and set 

X R {t) = xf R{T ^ k {t - T k ) + A$j(t), T k < t < T k+1 , k > 0, 

x R (T s - { X R (T k -) + ALi(T k ), if X? (T fc -) + ALi(T k ) >b, 

1 [ k) ~\b, ifxf(r fc -) + AL. 4 (r fc )<6, k - l > 

for all 1 < i < N where A^T^j(i), t>T k , is the local time at b of the process 

X R (t), t>T k . We remark that this is precisely the construction of X with 
the regulation and reflection at the barrier being added at the appropri- 
ate random times. Due to the definition of the processes X C)X)l and of the 
involved local time processes, X R and the corresponding process of regu- 
lations solve the regulated system (1.1), (1-2). The uniqueness of the weak 
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solution to (2.1), (2.2) implies that the weak solution of the regulated system 
is unique. 

It follows immediately that the gap processes Z(t), Z R {t) are well defined 
and unique in law. Their dynamics are given in the next lemma. 

Lemma 2.1. (a) The components of the gap process Z in the unregulated 
particle system satisfy 

dZi{t) = (6 i+ i -5 l )dt + a i+1 dp i+1 (t) - a x dpi(t) + dX i+1 (t) - dAi(i) 

i+i,i+2) (t) + dA( lj2 ) (i)) 
+ (Fi(Z(t-),AX(t)) - (Zi(t-) + A(X i+1 (t) - Ai(t)))). 

(b) The components of the gap process Z R in the regulated particle system 
are governed by 

dZ R {t) = 5 i dt + a l dp i {t) + dX i {t) + £(dAg_ 1)0 (t) - dAf i i+1) (t))l m} 
+ (dAg il) (t)-idAf 1)2) (i))l {i=1} 

+ (F R (Z*(t-),AA(t)) - (Z?(t-) + AA,(t))). 

Hereby, /3\(t), . . . ,f3^{t) are i.i.d. standard Brownian motions, X{t) = (Ai(t), . . . , 
Ajv(t)) is a Levy process whose components are i.i.d. pure jump Levy pro- 
cesses with the same distribution as L\(t) and independent of j3\{t),. . . ,f3^{t), 
A(j j +1 )(t) and A R i+1 ^(t) are the local times at of the processes Yi + \(t) — 

Yi(t) and Y^^t) - Y/ 1 ^), respectively, with A (A r iJV+1 ) (t) = A^ |iV+1) (i) = 0, 
and A^ ^(t) is the local time of the process Y R (t) at b. Finally, F and F R 
describe the value of the gap process after a jump as a function of its value 
before the jump and the jump of X . More explicitly, 

Fi(z,v)= min max(z jl _i+7 ?il ,...,z ii+1 _ 1 +7 ?ii+1 ) 
i<n<---<ji+i< N 

- min (zj-i + rjj) 

for all l<i<N -I, z = (zx, . . . , z N -i) E and rj = (r/i, . . . ,ri N ) & R N 

with the convention zq = and 

F R (z,i]) = min max(max(z, 1 + 77^,0), . . . ,m&x(zj i + %,0)) 
i<ji<-<ji<A r 

for all l<i<N, z = (z 1 ,...,z N ) € and n = (rji, . . . ,r] N ) G R N . More- 
over, the state spaces the processes are given by 

W N ~ 1 = {zi,..., z N -i\0 <Z!<z 2 <---< zjv-i} C 

and W N , respectively. 
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Remark. The contribution of the terms 

(Fi(Z(t-),AX(t)) - (Zi(t-) + A(X l+1 (t) - Ai (*)))), 1 < i < N - 1, 
and 

(F t R (Z R (t-),AX(t)) - (Z R (t-) + AAi(t))), 1 < i < iV, 

in the respective dynamics can be understood as follows. The contribution 
is nonzero if and only if one of the particles jumps and this jump changes 
the ranks of the particles. If this jump does not change the left-most particle 
in the unregulated system or does not involve a regulation in the regulated 
system, then these terms correspond to consecutive normal reflections of 
the process of gaps at faces of W N_1 or W N , respectively. More precisely, 
the gap process is normally reflected at the faces \zj x — Zj^i}, . . . , {zj 2 ~ l — 
Zj 2 }, if the inequalities Zj 1 < Zj 1+ i, . . . , Zj 1 < Zj 2 are violated by the jump 
and at the faces {zj 2 = Zj 2 -±}, . . . ,{zj 1+ i = Zj x }, if the inequalities Zj 1 < 
Zj 2 , . . . , Zj 2 -i < Zj 2 are violated by the jump. If the jump changes the left- 
most particle in the unregulated system, then the particles are relabeled and 
this term gives the change of the gaps due to relabeling. If a particle jumps 
below the barrier in the regulated evolution, then the particle configuration 
is regulated and the particles are relabeled. The term then gives the change 
of gaps due to both these operations. 

Proof of Lemma 2.1. (a) For any fixed t>0 let m : {1, . . . , N} — > 
{1, . . . , N} be a bijection such that 

(2.7) ^ ^irrW*") " ' ' ' " X ^\N)(t-)- 

In Section 3 of [3] the authors show that in the absence of jumps, that is, 
when L(t) = 0, it holds 

N 1 

(2.8) dYi{t) =J2 1 {* t {j)=i}dX j (t) + -(dA (i _ M) (t) - dA (M+1) (t)), 

i=i 

where Ar i i+ i)(t) are defined as in the statement of the lemma for 1 < i < N 
and we have set Ar 01 \(t) = 0. Moreover, equation (3.3) of [3] states that 
there exist i.i.d. standard Brownian motions /3i(i), . . . , Pn{P) such that 

(2.9) dYi (t) = 5idt + oi dPi (t) + \ (dA^i.i) it) - dA (iii+ i) (t)) • 

This yields immediately the claim of part (a) of the lemma if L(t) = 0. In 
the presence of jumps we define the pure jump processes Xi(t), . . . , Ajv(t) by 
AAj(i) = AL^-i^Jt). The latter are i.i.d. pure jump Levy processes with 
the same law as L\(t). Indeed, the jump times of A(t) and L(t) coincide 
and the law of AA(t) conditional on AX(t) ^ is the same as the law of 
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AL(t) conditional on AL(t) 7^ 0. Next, we recall that Tk, k > were defined 
by To = and as the jump times of the process L(t) for k > 1. Since X(t) 
coincides with the corresponding process in the absence of jumps for T^_i < 
t <Tk and any k > 1 by its construction, it remains to verify that AZj(Tfc) 
coincides with the jump given by the right-hand side of the equation in part 
(a) of the lemma for all 1 < i < N — 1 and k > 1. But this follows directly 
from the definition of F. 

(b) As in (a), we first treat the case Lit) = 0. In this case we will show 
that the processes Z R (t) = Y^(t) -Y^t), l<i<N, with Y R (t) = b follow 
dynamics corresponding to the dynamics in part (b) of the lemma. To this 
end, we set 

(2.10) Nj(t) = \{l<i< N\X R (t) = Xg)(t)}| 

for all 1 < j < N, t > 0. Next, we observe that Theorem 2.3 of [4] applied to 
the continuous semi-martingales X^t), . . . ,X^(t) yields for 1 < j < N, 

N j-1 
dY 3 R (t) = ^(N,(t)r\ x n )(t)=xm dX t R (t) + (Nj^r^dAfk^t) 

i=l 3 k=l 

-(NM' 1 E dAf m (t), 
k=j+i 

where j 2 \{t) is the local time of Y R (t) — Y R (t) at zero for 1 < ji < ji < 
N. Plugging (1.1) and (1.2) into the latter equation and applying the strong 
Markov property of (Z R (t), . . . , Z R (t)) to the entrance times of the set 

d £ R^ = {ze R^\dist(z,dM^) > e} 

for a fixed e > 0, one shows that Z R (t) = (Z R (t), Z R (t)) evolves as a 
Brownian motion 

(a 1 p 1 (t),a 2 P2{t) ~ o"i/3i(i), . . .,a N f3 N (t) - ctjv-i£jV-i(£)) 

with constant drift vector (61,62 — 5\, . . . ,5n — 6n-i) between an entrance 
time of the set d £ M.+ and the first hitting time of after that. This is due 
to an application of Knight's theorem in the form of [21], page 183, to the 
martingale parts of Y R (t), 1 < j < N. Consequently, letting e tend to zero, 

we observe that Z R (t) is a semi-martingale reflected Brownian motion in 
the orthant in the sense of [20, 25]. Now, Lemma 2.1 of [25] shows that 
the Lebesgue measure of the set {t > 0\Z R (t) £ cM+} is zero almost surely 
and Theorem 1 of [20] implies that the times for which either Nj(t) > 3 or 
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Nj(t) = 2 and Z\{t) = do not contribute to the dynamics of Y R {t) for all 
1 < J ' < N. Hence, the dynamics simplifies to 

dY R {t) = 5j dt + aj dfy (t) + {dAf 01) (t) - \ dAf 12) (t)) ■ l {j=1} 

+ |(dAg_ 1| . ) (i)-dAg J . +1) (t))-l { ^ 1} 

for all 1 < j < N. This yields immediately the statement of part (b) of the 
lemma in the absence of jumps. The general case follows by defining X(t) in 
the same way as in the proof of part (a) of the lemma and by making the 
next two observations. First, the dynamics of the process Z R {t) between the 
jump times of the process L(t) coincides with the dynamics of the system 
in the absence of jumps as a consequence of the construction of the process 
X R (t). Second, the jumps of the process Z R (t) coincide with the jumps of 
the right-hand side of the equation in part (b) of the lemma due to the 
definition of F R . □ 

2.2. Invariant distributions and convergence. In this section we inves- 
tigate the existence and uniqueness of invariant distributions of the gap 
processes in the two finite particle systems, as well as the convergence of the 
processes of gaps to the respective invariant distributions. We start with the 
proof of Theorem 1.2. 

Proof of Theorem 1.2. (a) (1) We first prove that the family Z(t), 
t > 0, is tight for any initial value z in W N ~ 1 . To do this, it suffices to show 
that the family Zi(t) = Zi(t) - Zi_i(i), t > 0, is tight for all 1 < i < N - 1 
where we have set Zo(t) = 0. To this end, we fix a 1 < i < N — 1, set y = Zi(0) 
and define the process D- on the same probability space as X by 

D y t (t) = y + (S i+1 - St)t + a i+1 /3 i+1 {t) - a^t) 

< 2 - n > 

+ XE lAA^-WI+A^tt), 

j=l 0< s <t 

where A D ^ (t) is the local time at of the process D^it) and /3\(t), . . . , /Sjv(i) 
and Ai(t), . . . , Ajv(i) are the same as in the dynamics of Z(t) given in Lem- 
ma 2.1(a). We note that AZ^t) < A-DV(t) for all t > and that after each 
time t > with Zi(t) = Df(t) the processes Zj and Df evolve in the same 
way until either the (i — l)st ranked particle and the ith ranked particle 
collide, or the (i + l)st ranked particle and the (i + 2)nd ranked particle 
collide, or there is a jump of (Xi(t), . . . , Xjy(t)). In particular, for any t > 
the accumulated local time at on the set {0 < s < t\Zi(s) = D\ (s) = 0} is 
the same for the processes and Putting these observations together, 
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we conclude that Zi(t) < Djit) almost surely for all t > 0. Hence, to prove 
that the family Zi(t), t>0 is tight, it suffices to show that the family D^{t), 
t > is tight. 

(2) Next, we fix an e > and find a C = C(e) > such that 



(2.12) P( sup D$(u)<C) >0. 



0<u<e 

We claim that the process Df(ne), n G N, is a recurrent Harris chain on 
E_l_ with respect to the set [0, C] in the sense of Section 5.6 of [10]. Indeed, 
the Harris property follows immediately from the Harris property of the 
corresponding process in the absence of jumps. Moreover, the law of large 
numbers for Levy processes and condition (1.8) show that almost surely 




(<S»+i - 5i)t + a i+ if3 i+1 (t) 



N \ 

*Mt) + Y, E l AA ;( fi )l 

j=l0<s<t J 



N 

L 0<s<l 



Hence, there exist stopping times < t\ < T2 < • ■ • tending to infinity with 
^ii T k) = for all k € N with probability 1. This follows by setting 

n =inf{n> l\DV(u) = 0}, 

r k+1 = inf{n > r k + l\D?(u) =0}, k > 1, 

and coupling Df(t), t> s with the Levy process 

£>V(s) + (S i+1 - 6i)(t -s) + cr i+1 (p i+1 (t) - A+i(s)) - <Ti(pi(t) - Pi(s)) 

N 

j=l s<u<t 

to conclude 

F({DUs) = 0}U{3t>s\Dy(t) = 0}) = l 

for all s > 0. In addition, we let Tfc(e) be the integer multiple of e which is 
closest to Tfc from above for all fcGN and observe that the strong Markov 
property of Di applied to the stopping times T k , k G N, (2.12) and the second 
Borel-Cantelli lemma imply that DNt^e)) < C for infinitely many k € N 
almost surely. This shows the recurrence of D^{ne), n G N. Moreover, by 
the results of Section 5.6c of [10], the chain is aperiodic and converges in 
total variation to its unique invariant distribution which we denote by v e . 
Furthermore, the uniqueness of invariant distributions of recurrent Harris 
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chains implies that u\ = V\j2 = = • • • • Next, we fix a ( > and at«>0 
and define the function f w :R+ — y M + by f w (v) = e~ w ' v . Since Di is a Feller 
process [this can be shown along the lines of step (3) below], its semi-group 
of transition operators is strongly continuous on the space of continuous 
functions on R_|_ vanishing at infinity (see, e.g., Theorem 17.6 in [15]). Thus, 
denoting the semi-group of transition operators corresponding to Di(t), t > 
0, by P A (t), t>0, we can find aQsN such that 

(2.13) V0 < t < 2~ Q : \\P Di (t)f w - /JU < C- 

Moreover, for each t > we let n(t) be the largest integer such that n(t) ■ 
2~® < t. All in all, setting /j, Vi t = 5 y P Di (t) and writing //(/) for f fd/j,, we 
get for any t > 0, 

\f*y,t(U) ~ v l(fw)\ 

< \Vy,t(fw) ~ Vy,n(t)2-Q(fu,)\ + \Vy,n{t)2-Q (fw) ~ Vl(fw)\ 

= \Vy,n{t)2-Q{P Dl {t-n(t)2- Q )f w - f w )\ + \/J.y jn ( t )2-Q(fw)-Mfw)\ 

< C + \f*y,n(t)2-Q(fv>) ~ u l{fw)\- 

Taking first the limit t — > 00 and then the limit £ J, 0, we conclude that 
Vy,t(fw) converges to fi(f w ). Since w > was arbitrary, it follows that the 
Laplace transforms of the measures 5 y P Di (t) converge point-wise to the 
Laplace transform of v±, hence, the measures 5 y P Di (t) converge weakly to 
v\. Thus, the family DV(t), t>0, converges in law to v\ and is, therefore, 
tight as claimed. 

(3) Next, denote by P(t), t>0, the Markov semi-group of operators cor- 
responding to the Markov process Z. The tightness of Z(t), t > implies 
that the family of probability measures ~ J * 5 z P{s) ds, t > 0, is tight. Hence, 
we can find a sequence t\ < t2 < ■ ■ ■ tending to infinity such that the weak 
limit 

(2.14) u= lim — f n (5 z P(s))ds 

n->oo t n J 

exists and is a probability measure on W N ~ l . We claim that v is an invariant 
distribution of Z. To this end, we first remark that Z is a Feller process. For 
the evolution without jumps this is a consequence of the Feller property of 
the reflected Brownian motion (see Theorem 1.1 and the following remark 
in [26]) and the observation that the process (Zi(t), . . . , Zjv-i(t)) can be 
viewed as a reflected Brownian motion in the sense of [26] (see [14] , Section 
5). In our case the Feller property can be seen as follows. Let z° be an 
arbitrary point in W N ~ l and {z n )'^ =l be a sequence in W N_1 converging to 

71 

it. Moreover, let Z z and Z z , n > 1, be the gap processes with initial values 
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z° and z n , n> 1, respectively. We need to show that for any fixed t > the 
random vectors Z z (t) converge in law to Z z (t). To this end, we note that 
the gap process Z can be constructed by first generating the sequence of its 
jump times T\ < T2 < ■ ■ ■ and the corresponding jumps AA(Ti), AA(T2), . . . 
of X(t) = (Ai(i), . . . , Ajv(i)) an d then defining Z conditional on these choices 
by 

00 

(2.15) Z(t) = £ l {Tk < t<Tk+l} Z^ziT k - )A HT k)) , k{t _ T ^ 

k=0 

where we have set To = 0, Z(0—) = Z(0), F is the continuous function as in 
part (a) of the Lemma 2.1 and Z c ' F ( z ( Tk ~^ AX ( Tk ^' k , k > 0, are independent 
gap processes of unregulated evolutions without jumps with respective ini- 
tial values F(Z(T)~— ), AA(T^)), k > 0. Due to the Dominated Convergence 
theorem it suffices to show that the law of Z z (t) conditional on a real- 
ization of A converges weakly to the law of Z z (t) conditional on the same 
realization of A. But in view of the representation (2.15), this can be shown 
by using induction over the unique value of k for which t € [T^, Tfc+i) and 
the Feller property of Z c '' ,k . 

Moreover, for any t > we have 

1 ft n 1 ftn+t 

(2.16) lim — / (S z P(s))ds= lim — / (6 z P{s))ds, 

n->oo t n J n->co t n J t 

where the limits are taken in the weak sense. This is a consequence of the 
fact that the total variation norm of 4- J*q™ (8 z P(s)) ds — j- //™ +< (<5 2 -P(s)) ds 
is bounded above by |^ for all n with t n > t. Hence, by the Feller property 
of Z 

[ fd(uP(t))= lim / P{t)fd(^f n {5 z P(s))ds 
= lim I fd(±- [ tn+ \5 z P(s))ds 



J 

Jw 



fdv 

for all continuous bounded functions /. Thus, v is an invariant distribution 
of Z. 

(4) We now prove that v is the only invariant distribution. To this end, 
consider the process Z(n), n € N. We claim that it is a recurrent Harris chain 
in the sense of Section 5.6 in [10] . Indeed, in the absence of jumps the set {t > 
0\Z(t) G dW N ~ 1 } has Lebesgue measure zero almost surely and the process 
Z evolves as a Brownian motion (o^/^W — &iPi{t), ■ ■ ■ , crN^N(t) — ci/3i(i)) 
with constant drift vector (82 — . . . ,5n — 5i) in the interior of W 1 ^^ 1 as we 
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have seen in the proof of Lemma 2.1(a). Since the covariance matrix of the 
latter Brownian motion is nondegenerate and there is a positive probability 
that Z has no jumps in the time interval [0,1], we conclude that Z(n), 
n £ N, is a Harris chain on W N ~ l . Noting that v is an invariant distribution 
for this chain, we conclude that the chain must be recurrent. Thus, it has a 
unique invariant distribution (see [10], Section 5.6). Since any other invariant 
distribution of Z is an invariant distribution of Z(n), n € N, it has to coincide 
with v. 

(b) (1) Part (b) of the theorem can be established by using the technique 
of the proof of part (a). However, we present here a softer argument based 
on a monotonicity property special to the regulated system. The latter was 
inspired by the proof of Lemma 1 on page 162 of [7] which deals with a similar 
discrete time problem. The crucial idea is to introduce a family of processes 
indexed by a > such that each of them evolves as the regulated system 
after time —a and for any two indices a±, a<i the corresponding processes 
are driven by the same Brownian motions and pure jump Levy processes 
after time max(— a\, —02)- To construct such a family we start by defining 
auxilliary independent standard Brownian motions /3 mj i(t), . . . , /3 m ,jv(t) and 
independent i.i.d. pure jump Levy processes A m) i(i), . . . , X m ^(t) of the same 
law as Ai(t) for all m £ Z_ on an extension of the probability space (O, J 7 , P) 
on which X R was defined. Moreover, we let m(t) be the largest integer which 
is less than or equal to t and l(t) = t — m(t) for all tel. Finally, we can use 
the notation to define a family of processes V a on W N indexed by a > 
with the desired properties by setting V a (t) = b ■ 1 = (6, . . . , b) for t < —a, 

dV?{t) = Si dt + <Ti dp mm (l(t)) + d\ m (t),Mt)) 
+ UdA? i _ 1>i) (t)-dAl i+1) (t))-l {im 
+ (dAf 0il) (t)-idAf 1)2) (i))-l {i=1} 

+ (F?(y a (t-), AA m(t) (Z(t))) - (V?(t-) + AX mm (l(t)))) 
for —a <t<0 and 

dvf(t) = s l dt + dfc(t) + d\i(t) + U dA (i-i,i) (*) - d %+D (*)) • Him 

+ (dAf 0)1) (t)-idAf li2) (t)).l {i=1} 

+ (F l R (V a (t-),AX(t)) - (V?(t-) + AAj(t))) 

for t > 0. Hereby, notation is as in Lemma 2.1(b) with the local times defined 
with respect to V a instead of Y . We note that the processes V a are defined 
in such a way that the law of V a (t) for t > —a is the law of the ordered 
particle configuration in the regulated evolution at time t + a, started in 
6-1. 
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(2) Moreover, we have for V ai , V a2 with a\ < a<i the inequalities 

(2.17) V? (t) < (t), . . . , (t) < (t) 

for all t > —a.\. Indeed, this is clear for t = —a±. Furthermore, the inequal- 
ities are preserved under the jumps of the two processes because the jump 
times and the jump sizes before the relabeling and the regulation of the par- 
ticle configuration are the same for both processes. In addition, started with 
each stopping time at which a nonempty set J C {1, . .. ,N} of coordinate 
processes of V ai and V° 2 coincide, these coordinate processes evolve in the 
same way until one of the particles with rank in J collides with a particle 
with rank not in J or until there is a jump which changes the rank of at least 
one particle whose rank was originally in J. By distinguishing the two cases, 
one concludes that the regulated evolution preserves the component-wise 
order < on configurations in W N . Hence, for all 1 < i < N we may define 

(2.18) V°°(t)=t lim V?(t) 

a— >oo 

as an element of [b, oo]. 

(3) Next, we claim that the random vector y°°(0) is finite almost surely. 
To this end, we observe that the process Z R (t), t > 0, started at Z R (0) = 
converges in distribution to V°°(0) — b ■ 1 for t — > oo, because 

(2.19) Z R (t)±(V 1 t (0)-b,...,Vh(0)-b) 

and the right-hand side converges to V°°(0) — b-1 almost surely. In addition, 
we have for all 1 < i < N and all n 6 N that 

(2.20) Zf\n) = £ Zf(n) < £ Df{n) 

3=1 3=1 

with Zf(t) = Zf(t) - Zf_ x (t) for all t > 0, 1 < j < N, Z§{t) = 0, 

N 

Df t (t) = 6 1 t + a 1 /3 1 {t)+J2 |AA fc ( S )|+A D ?(t), 



k=10<s<t 

N 



D R (t) = (6 j -S j - 1 )t + a j p j (t)-a j - 1 f3^ 1 (t) + J2 E \ AX ^)\ + A^(t), 

fc=10<s<i 

2 < j < N and A D ?(t) being the local time of Df{t) at zero for l<j<N. 
We now use condition (1.9) and follow the lines of step (2) of the proof of 
part (a) to conclude that D R {t) converges in distribution to an almost surely 
finite random variable in the limit t — > oo for each 1 < j < N. Thus, for all 
l<i<N, the sequence Y^j=i^f ( n )' n £ N, is tight. Hence, the sequence 
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Z R (n), n £ N, is also tight, so the limit in distribution of Z R (n), n € N, 
which is V^°°(0) — 6, must be almost surely finite for all 1 < i < N. 

(4) Next, we prove that Z R is a Feller process. To this end, it suffices to 
prove the Feller property in the absence of jumps and to follow the lines 
of step (3) in the proof of part (a) to extend it to the general case. But 
in the absence of jumps, the process Z R is a reflected Brownian motion in 
W N in the sense of [26] as we have shown in the proof of part (b) of the 
Lemma 2.1. Thus, by Theorem 1.1 and the following remark in [26], the 
process Z R is a Feller process in the absence of jumps. The Feller property 
and the convergence in law of Z R (t), t > to V°°(0) — b ■ 1 imply that the 
law of y°°(0) — b • 1 is an invariant distribution of the process Z R . The 
uniqueness of the invariant distribution can be deduced from the uniqueness 
of the invariant distribution of the chain Z R (n), nSNby following step (4) 
of the proof of part (a). □ 

Now, we are able to prove Theorem 1.3 which deals with the convergence 
of the gap processes to the respective invariant distributions. 

Proof of Theorem 1.3. In the course of the proof of Lemma 2.1 we 
have seen that in the absence of jumps, the processes Z and Z R are reflected 
Brownian motions in the sense of [26]. Moreover, as observed in the proof 
of Theorem 1.2, their respective covariance matrices are nondegenerate. In 
the presence of jumps this shows the Harris property of the processes Z(t), 
t > 0, and Z R (t), t > 0, as well as the Harris property and the irreducibility 
of the chains Z(n), n 6 N, and Z R (n), n E N. The existence of invariant dis- 
tributions (Theorem 1.2) implies that the processes Z(t), t>0, and Z R (t), 
t>0, are positive Harris recurrent in the sense of [18]. Thus, Theorem 6.1 
of [18] shows that Z(t), t > 0, and Z R (t),t > 0, converge in total variation 
to their respective invariant distributions. □ 

3. Systems of infinitely many particles. 

3.1. Existence and uniqueness in law of the processes. From now on we 
let / = N. Due to the fact that the processes Li(t),L2(t), . . . are i.i.d., an ap- 
plication of the second Borel-Cantelli lemma shows immediately that in both 
evolutions infinitely many jumps of the particles will occur almost surely on 
each nonempty time interval, if L\{t) is not identically equal to zero. Hence, 
the proof of existence of the finite particle systems does not carry over to 
the infinite case. Also, the proof in [19] for the case of the unregulated sys- 
tem without jumps, which uses Girsanov's theorem for Brownian motion, 
cannot be applied here, since it would prove the existence of the solution to 
the unregulated version of (1.1) in the case of a certain nontrivial depen- 
dence structure between the processes B\(t), B2(t), . . . and L±(t),L2(t), 
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Instead, we prove the existence and the uniqueness in law of the infinite 
particle systems by a bound on the tail of the distribution of the running 
supremum of an integrable Levy process. 

Proposition 3.1. If the initial configuration Xi(0), X 2 (0), . . . of the 
particles satisfies 



oo 




for all t > and y G M ; then the unregulated version of (1.1) has a unique 
weak solution and the corresponding ordered particle system Yi(t),Y2(t), . . . 
is well defined for all £ > 0. 

In particular, condition (3.1) is satisfied if there are constants 71 > 0, 
72 G R with 

(3.2) Xi(0)>7ii + 72 , i>l. 

The same statements are true for the regulated system (1.1), (1-2). 

Proof. (1) We prove the proposition only for the unregulated system, 
since the assertion for the regulated system can be shown in the same way by 
setting M = 1. To this end, we assume (3.1) and introduce the probability 
space 

( 11 x n B ' L , jfW*) <g> f b > l , (g) pW*) ®f B A. 

\n,x) (N,x) (N,x) ' 

Hereby, the products are over {(N, x)\N > l,x € R , x\ < ■ ■ ■ < x^}, on each 
(Jl^), J<^»*),pW*)) the unregulated system X N > Xl >-> XN with N particles, 
initial configuration x\ < X2 < ■ • ■ < x^ and parameters S\, . . . , 5n, o"i, . . . , g\ 
are defined and (Q B,L , J rB,L ,P B,L ) is a probability space on which the i.i.d. 
Brownian motions Bi(t),B2(t), . . . and the i.i.d. pure jump Levy processes 
Li(t) , L2(t) , . . . independent of the Brownian motions are defined. We call 
the product space above (f2°°, J-°°, P°°). Next, we recall the definition of the 
constant M introduced in Assumption 1.1, set Gi(t) = a\Bi{t) + Li(t) for 
all t > and i & I for the sake of shorter notation and define the process X, 
the sets Aq, C A\ C • • • and the stopping times = po < p\ < ■ • ■ inductively 
by 

A k = {% > 1|31 < j < M, < s < p k : Xi(s) = X (j) (s)}, 

X i (s) = X i (s-pk), iEA k ,p k <s<p k+1 , 

Xi(s) = Xi(p k ) + 8 M (s - p k ) + Gi(s) - Gi(p k ), i(£ A k ,p k <s < p k+1 , 
p k+ i = inf{s > p k \31 <j< M,i<£ A k :Xi(s) = X {j) (s)}. 
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Fixing a t > we observe that by condition (3.1) and the first Borel-Cantelli 
lemma there are almost surely only finitely many particles which visit the 
interval (— oo,y\ up to time t for any fixed y £ M. This follows by setting 
y = y + maxj>i|<5j| • t in (3.1) and recalling that the set of collision times 
of distinct particles in the finite unregulated particle system has Lebesgue 
measure zero almost surely (see the proof of Lemma 2.1). Furthermore, we 
note that by choosing a large enough y we can make 
M 

Vp(x,-(0) + sup (ai J Bi(s) + L 1 (s))+max|5 i |-t>y N ) 

V 0<s<t *>1 ' 

j=l 

as small as we want. We also note that in the regulated case this expression 
should be replaced by 

sup (Xi (0) + Sis + <txB x (b) + Lx(s) + A (1 6) (s))>y 

0<s<t 

This observation implies that there exists a K = K(ca) such that K < oo, 
Pk < t < PK+x and the sets . . . , Ak are finite almost surely. Thus, X is 
well defined on [0,t] for almost every w € Since i > was arbitrary, we 
have shown that X is well defined on [0, oo) for almost every uj 6 By its 
construction, X is a weak solution to the unregulated version of the system 
(1.1). Moreover, at any time t > the ordered particle system Yi(t),Y2(t), . . . 
is well defined, since there are finitely many particles on each interval of the 
form (— oo,?/] almost surely [due to (3.1) with y = y + maxj>i|5j| • t and the 
first Borel-Cantelli lemma]. 

(2) For the uniqueness part let X' be a different weak solution to the 
unregulated version of (1.1). Then we can define inductively the sets A' Q C 
A[ C • • • and the stopping times = p' < p' ± < ■ ■ ■ by 

(3.3) A' k = {i > 1|31 < j < M, 0<s<p' k : X[(s) = X[ j} (s)}, 

(3.4) p' k+1 = mf{s > p' k \31 < j < M, i $ A k : Xfa) = X[ j} (s)}. 

Due to the uniqueness of the weak solution to the unregulated system (1.1) in 
the case of finitely many particles, the joint distribution of p' , p±, . . . and X' 
on [p' Q , p'-j], [p'i,P2], ■ ■ ■ has to coincide with the joint distribution of po,Pi, ■ ■ ■ 

and X on [po,/?i] ; [pi ; P2]; Thus, the law of X' is the same as the law of 

X. 

(3) To prove that (3.2) implies (3.1) we first observe that 

(3.5) E[max(|<7iBi(t) + Li(t)|,l)] < oo 

for any fixed t > 0. Noting that max(|x|,l) is a nonnegative continuous 
submultiplicative function (see Proposition 25.4 of [23]), we conclude from 
the Theorem 25.18 of [23] that 



(3.6) E sup |<7iBi(a) + Li(a) 

H)<s<t 



< oo 
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for any fixed t > 0. Thus, (3.2) implies (3.1). □ 

3.2. Tightness of the infinite regulated system. We can now prove Propo- 
sition 1.4 which guarantees the tightness of the gap process in the infinite 
regulated evolution. 

Proof of Proposition 1.4. We start with the observation that for 
each i £ / the process X^it), t > and the corresponding process of regu- 
lations Ri(t), t>0, solve the system of stochastic differential equations 

dX^(t) = 5\ dt + eri dBi(t) + dLi(t) + dR^t), 

dRi(t) = (b- X?(t-) - AL i (i))l {xf(t _ )+Aii(t)<b} + dA (i>6) (i). 

This is due to the fact that the drift and the volatility sequences are con- 
stant and that the set of collision times of distinct particles has Lebesgue 
measure zero almost surely. The latter statement is a consequence of the 
corresponding property of the regulated system with finitely many particles 
[see the proof of part (b) of Lemma 2.1] and the construction of the infinite 
particle systems (see the proof of Proposition 3.1). Next, for each % £ / we 
introduce the process 

(3.7) H^(t) = Xf (0) + 5 1 t + a 1 B l (t)+ £ |AL,( S )| + A H *(i), 

0<s<t 

where A^(i) is the local time of H-* at b. Using condition (1.10) and arguing 
as in the steps (1) and (2) of the proof of Theorem 1.2(a), we conclude that 
X^(t) < H,f (t) for all t > and i £ I almost surely and, in addition, that 
for each i£ I the process H^(t), t>0, converges in law to an almost surely 
finite random variable H^(oo) whose law does not depend on X^(0) and i. 
Moreover, due to the independence of the processes H^, i€l, the random 
vector (H^(t),...,Hf-(t)) converges m distribution to (tff (oo), . . . , Hf (oo)) 
for any j € I where H^(oo), i > 2, are chosen as independent copies of 
Hi(oo). Finally, the chain of inequalities 

(3.8) < Zf (f) < max Xf (t) -b< max #f (i) - b 

for all t > shows that the family Z^{t), t > 0, is tight for all j £ /. This 
yields the tightness of Z R (t), t > 0, on with the product topology. □ 
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